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We study electronic structure of hole- and electron-doped Mott insulators in the two-dimensional Hubbard 
model to reach a unified picture for the normal state of cuprate high-Tc superconductors. By using a cluster 
extension of the dynamical mean-field theory, we demonstrate that structure of coexisting zeros and poles of 
the single-particle Green's function holds the key to understand Mott physics in the underdoped region. We 
show evidence for the emergence of non-Fermi-liquid phase caused by the topological quantum phase transition 
of Fermi surface by analyzing low-energy charge dynamics. The spectra calculated in a wide range of energy 
and momentum reproduce various anomalous properties observed in experiments for the high-Tc cuprates. Our 
results reveal that the pseudogap in hole-doped cuprates has a d-wave-like structure only below the Fermi level, 
while it retains non-d-wave structure with a fully opened gap above the Fermi energy even in the nodal direction 
due to a zero surface extending over the entire Brillouin zone. In addition to the non-d-wave pseudogap, the 
present comprehensive identifications of the spectral asymmetry as to the Fermi energy, the Fermi arc, and the 
back-bending behavior of the dispersion, waterfall, and low-energy kink, in agreement with the experimental 
anomalies of the cuprates, do not support that these originate from (the precursors of) symmetry breakings such 
as the preformed pairing and the d-density wave fluctuations, but support that they are direct consequences of 
the proximity to the Mott insulator. Several possible experiments are further proposed to prove or disprove our 
zero mechanism. 

PACS numbers: 71.10.Hf; 74.72.Kf; 79.60.-i 



I. INTRODUCTION 

Anomalous behaviors of high-Tc cuprates observed in the 
normal metallic state above hold the key not only to under- 
standing the mechanism of the superconductivity but also to 
a possible manifestation of an unexplored metallic phase dis- 
tinguished from the Fermi liquid.' Extraordinary electronic 
structure is induced by a small density of carrier doping into 
the Mott insulator. Angle-resolved photoemission spectro- 
scopies (ARPES) have in fact revealed detailed anomalies of 
the normal-state spectra, such as momentum-dependent exci- 
tation gap (pseudogap), a truncated Fermi surface (Fermi arc), 
and kinks in the dispersion.^ 

Toward the understanding of the anomalous metals, es- 
pecially the pseudogap formation, many theoretical propos- 
als have been made so far.^ The proposals include a Cooper 
paring without phase coherence,"* and hidden orders or its 
fluctuations competing with the superconductivity, such as 
antiferromagnetism,^^ charge or stripe orderings,^ and d- 
density wave.^ '" Mechanisms attributing the origin of the 
pseudogap to a direct consequence of the proximity to the 
Mott insulator have also been proposed.""'^ 

Among the theoretical efforts, recent development of 
the dynamical mean-field theory (DMFT)'* and its cluster 
extensions'^ '^ has enabled studies on dynamics of micro- 
scopic models without any ad hoc approximation. In par- 
ticular, studies on the two-dimensional (2D) Hubbard model 
using the cluster DMFT (CDMFT) have offered many use- 
ful insights into the electronic structure of cuprates, by iden- 
tifying the pseudogap, Fermi arc,'^''^''''"^^ and high-energy 
kink^^ in the calculated spectra. The CDMFT is specifically 
suited and powerful for this problem because of its nonper- 



turbative framework, namely, it is based on neither weak nor 
strong coupling expansions. Furthermore, it takes account 
of short-range spatial coiTelations within a cluster explicitly. 
These are big advantages in exploring momentum-resolved 
dynamics in the intermediate coupling region, which is rel- 
evant to physics of the cuprates. In particular, recent CDMFT 
studies on doped Mott insulators have revealed emergence of 
non-Fermi-liquid phases, characterized by unexpected coex- 
istence of two singularities at the Fermi level, one characteris- 
tic to the weak-coupling and the other to the strong-coupling 
regions. '3.'5.26.29.30 

The singularity characterizing weakly interacting metals, 
namely, Fermi liquids, is a pole of the single-particle Green's 




FIG. 1. (color online). Schematics of the sign of ReG for (a) normal 
metals and (b) correlation-driven insulating states, k is the momen- 
tum and o) is the energy measured from the Fermi level. 
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function G. Energy dependence of its locus in the momentum 
space determines band dispersions, particularly, the Fermi sur- 
face at the Fermi level. On the other hand, a crucial singularity 
in the strong coupling region is a zero of G. The situations are 
illustrated in Fig. 1, which schematically shows how the real 
part of G changes the sign in the energy-momentum space. 
Because ReG must be positive at high energy while negative 
at low energy, it has to change its sign at least once between 
these two regions. The sign change also has to take place in 
the Brillouin zone at the Fermi level, say, between the zone 
center and the zone boundary, if the dispersion far from the 
Fermi level ensures positive ReG in some part of momenta 
and negative ReG in the other part. In the normal metals, the 
sign change occurs at the poles of G (i.e., the band disper- 
sions), where ReG goes to — oo on the lower side of the band 
and comes back from -hoo on the upper side [Fig. 1(a)]. On 
the other hand, what happens when correlation effects induce 
a gap in the band at the Fermi level? Because of the absence 
of poles inside the gap by definition, there is no way for ReG 
to change its sign at the Fermi level except for getting through 
zeros [Fig. l(b)].^''^^ ReG = means the divergence of the 
self-energy, so that a zero of G is a singular point of the self- 
energy. 

The non-Fermi liquids discovered in CDMFT show a co- 
existence of these two characteristics — poles and zeros of 
G at the Fermi level. '■'■'^ In a lightly hole-doped region, the 
poles form a hole-pocket Fermi surface around the nodal di- 
rection [from k = (0, 0) to (tt, it)], while the zeros form a 
surface enclosing (vr, tt). It has been suggested that the Fermi 
arc emerges because of such peculiar electronic structure: The 
pocket loses the spectral intensity on the side closer to the zero 
surface, leaving the other side as an arc.'^''^'^'^ Similar mecha- 
nisms were proposed, based on an assumed functional form of 
Green's functions^^'^"* or a weakly-coupled chain model. ^^'^^ 
The pseudogap is also characterized by the same zero surface 
crossing the Fermi level. In the previous paper^^ we clar- 
ified the structure of poles and zeros in the whole energy- 
momentum space for hole-doped and undoped Mott insula- 
tors. We thus proposed a unified picture for understanding 
various puzzling features of Mott physics in the light of re- 
constructions of pole-zero structure under progressive doping 
from the Mott insulator to the Fermi liquid in the overdoped 
region. 

The aim of the present paper is to show how and to what ex- 
tent we can understand experimental findings within the 'pole- 
zero mechanism' . We implement CDMFT calculations at zero 
and finite temperatures, clarifying the relationship between 
zero-temperature pole-zero structure and the spectra observed 
in experiments at finite temperatures. We find that a number 
of non-Fermi-liquid aspects, i.e., spectral asymmetry as to the 
Fermi level, back-bending or incoherent feature of the disper- 
sion, high- and low-energy kinks, and Fermi arc or pockets in 
the hole- and electron-doped cuprates, are comprehensively 
understood within the pole-zero mechanism. In addition to 
a coherent picture for various experiments, our numerical re- 
sults predict a distinctive feature of the pseudogap: The nu- 
merical data support a fully opened gap above the Fermi level 
even in the nodal direction, which is incompatible with sce- 



narios on the basis of the rf-wave gap in the zero-temperature 
limit. Our result offers a mechanism distinct from the scenar- 
ios based on the preformed pairing or d-density-wave order. 

The paper is organized as follows. In Sec. II we intro- 
duce the 2D Hubbard model and present the central ideas 
of the CDMFT. We also describe the essence of numerical 
solvers for the effective cluster problem that we used in the 
present study. The numerical results are presented and dis- 
cussed in comparison with various experiments in Sec. III. 
From the low-energy pole-zero structure, we derive a sim- 
ple interpretation of several unusual features in hole-doped 
cuprates (Sec. IIIA-IIIC). We also find several anomalies in 
the band dispersion, which are consistent with ARPES data 
(Sec. IIID-IIIF). In Sec. IIIG we discuss the results of the 
Fermi surface in comparison with the recent ARPES obser- 
vation of the hole pocket. The comparison is further made 
in Sec. Ill H by changing the next-nearest-neighbor hopping 
to reach more quantitative understanding of the experimen- 
tal results. In Sec. IIIH, III I, and III J, we extend our the- 
ory to electron-doped cuprates and find good agreements with 
ARPES data. We summarize our results and make concluding 
remarks in Sec. IV. 



II. MODEL AND METHOD 

As a simplest model for high-Tc cuprates we take the Hub- 
bard Hamiltonian, 

H = ^ e(k)4^Ck(7 - ii^Uia + Ui-fTiii, (1) 

kcr ia i 

on a square lattice. Here Cko- (f^kcr) annihilates (creates) an 
electron of spin a with momentum k ~ [k^, ky), cia {c\^) is 
its Fourier component at site i, and = cl^Cia- U repre- 
sents the onsite Coulomb repulsion, ^, the chemical potential, 
and 

e(k) = — 2t(cos kx + cos ky) — At' cos kx cos ky, (2) 

where t (t') is the (next-)nearest-neighbor transfer integral. 

Based on a first principles calculation, the value of t was 
estimated to be ^ 0.4eV for La2Cu04.^^ —t'/t is con- 
sidered to be ~ 0.2 and ~ 0.4 for La2-xSrxCu04 and 
Bi2Sr2CaCu208+i, respectively. We adopt U = 8t or 12t, 
which are realistic values for the cuprates and indeed repro- 
duce the Mott insulating state for undoped case. 

In the CDMFT^^ we map the system (1) onto a model con- 
sisting of an TVc-site cluster C and bath degrees of freedom B. 
The bath is determined in a self-consistent way to provide an 
Nc X Nc dynamical mean-field matrix 5o(*w„) at temperature 
T(= 1/13), where cj„ = (2n + l)7rT. 

After the self-consistency loop converges, we calculate a 
quantity defined on the original lattice from those on the 
cluster, Q'^}^ This periodization procedure is based on the 
Fourier transformation truncated by the cluster size N^, 

Q''(k) = ^^[0%e*''---, (3) 
" ijec 
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where k is defined on the entire Brillouin zone of the original 
lattice and is the real-space vector connecting two cluster 
sites i and 

The truncation by a small TVc gives a good approximation to 
the thermodynamic limit of A'c — >^ oo if Q*^ is well localized 
within the cluster. In reality it is difficult to find a quantity 
which is short ranged in the entire parameter range of the in- 
teraction strength and the doping concentration. Therefore we 
need to choose an appropriate quantity to periodize, according 
to situations. For example, Q = E (E: self-energy) is a good 
choice for weakly interacting normal metals/^ but it becomes 
highly nonlocal and long ranged in the strong coupling region 
(e.g., in the Mott insulator). This is due to an appearance of 
zeros of G, i.e., poles of E in the momentum space (see AP- 
PENDIX A for further detail). On the other hand, in the Mott 
insulator Q = G is more appropriate because it is nearly local 
in the strong coupling regime.^^ 

Another choice for the periodization is the cumulant 



M = [ia;„-F/u-E]" 



(4) 



It was pointed out that the cumulant periodization works well 
in a wide range of U including the strong coupling, be- 
cause it is similar in the functional form to the atomic Green 
function.'^ Moreover, the periodization by the cumulant has 
an important feature, i.e., it can describe both poles and ze- 
ros of G at the same time, while the periodization by using E 
(G) describes only poles (zeros). This opens up the intriguing 
possibility of exploring the coexistence of poles and zeros at 
the Fermi level, as substantiation of anomalous metals. 

In the present study, we adopt the cumulant periodization 
scheme, Q = M, to investigate metals in the vicinity of the 
Mott insulator. In fact, M is highly local even in the doped 
metallic states. In APPENDIX B we demonstrate this local 
nature by CDMFT calculations for A^c = 4 x 4 cluster; we 
find that is nearly locaUzed already within the inner 2x2 
cluster for a parameter region relevant to the present study. 
After obtaining the lattice cumulant through Eq. (3), we 
calculate the self-energy E^, Green's function G^, and spec- 
tral function on the original lattice with 



GL(k, uj)=[oj + n- e(k) - EL(k, w)] , 



(5) 



A. Exact diagonalization method 

Although the pseudogap state in cuprates is experimentally 
detectable only above T^, it is still significant to elucidate the 
nature of the pseudogap state in the zero temperature limit, by 
assuming the paramagnetic metal for the doped Mott insula- 
tor This is a circumstance similar to the normal metal, where 
the concept of the Fermi liquid justified only in the zero tem- 
perature limit in the strict sense has proven to be fruitful. For 
this purpose of clarifying the zero temperature limit, we em- 
ploy the Lanczos ED method^^ for the cluster problem. 

In this scheme we take a large but finite value of pseudo 
inverse temperature /3' (= 100/t or 200/t throughout the pa- 
per), which practically represents the ground state with an en- 
ergy resolution corresponding to 1//3'. We then represent the 
dynamical mean field qq with a finite number A'b of bath de- 
grees of freedom, which constitute together with C sites the 
effective Hamiltonian to be diagonaUzed. We take A^b = 8 
throughout the paper 

The optimization of B sites is done by minimizing the dis- 
tance function defined by 



d. 



EE 



-w„/t 



(7) 



where 50. is the non-interacting Green's function for the ef- 
fective Hamiltonian and we have introduced the exponential 
weight factor e~'^"/* with w„ = (2n -|- l)7r//3' to reproduce 
more precisely the important low-frequency part. We have 
examined several other types of distance functions and con- 
firmed that qualitative feature of the results obtained in this 
paper does not depend on the choice. 

An advantage of the Lanczos method is that Green's func- 
tion is obtained as a function of real frequency w by the 
continued-fraction expansion. 
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and 



(8) 



A^{k,io) = JmG^{l!i,Lj). 



(6) 



In the following calculations, we employ an iVc = 2 x 2 clus- 
ter, and concentrate on the paramagnetic metallic solution. 

We numerically solve the effective cluster problem by 
means of the exact diagonalization (ED) method at T = and 
the continuous-time quantum Monte Carlo (CTQMC) method 
at T > 0, as we elaborate below. We hereafter omit the super- 
script L in E^, G^ and A^. 



where |0) is the ground-state vector and the coefficients af'^ 
and bf'^ are the elements of the tridiagonal matrix appearing 
in the Lanczos algorithm.^^ We take account of up to 2000th 
order in the expansion. A small positive rj is introduced to sat- 
isfy the causality. In principle, r] is taken to be infinitesimal, 
but in practice, it is useful to consider 77 as a parameter, which 
serves as a resolution in energy or mimics an infinite-size ef- 
fect not incorporated into the ED calculation. 

In the following study, we pursue a further benefit of the pa- 
rameter r]: We use as a mimic of the source of incoherence in 
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the electronic structure, such as thermal or impurity scattering. 
This is based on the observation that rj does not substantially 
change the location of poles and zeros of G but changes only 
the sharpness of these singularities. '^'^"^ Indeed, as long as rj is 
sufficiently smaller than the typical energy scale of the system, 
the location of the poles and zeros is virtually determined only 
by a^'^ and bf'^ in Eq. (8), which are calculated in the self- 
consistency loop performed on the Matsubara-frequency axis 
and independent from rj. Thus, 77 provides an opportunity to 
get insight into how the electronic structure at finite tempera- 
tures evolves from the pole-zero structure at zero temperature. 
The effect of 77 is confirmed by a direct comparison with the 
finite-temperature results obtained by QMC introduced below. 
Note that this smearing technique by rj is very useful because 
it is difficult for QMC to obtain the precise electronic struc- 
ture at real frequencies because it requires an analytic contin- 
uation of the numerical data. We use the smearing technique 
in Sec. Ill to compare the CDMFTh-ED results with ARPES 
ones. 



B. Quantum Monte Carlo method 

In order to discuss thermal effects directly, we implement 
QMC calculations for the effective cluster problem. Since the 
scheme takes into account infinite bath degrees of freedom, 
the results complement the limitation in the bath size in the 
CDMFTh-ED results as well. We adopt the algorithm based 
on a weak-coupling series expansion and auxiliary-field trans- 
formation. The idea was first proposed by Rombouts et al.f"^ 
applied to DMFT by Sakai et al.f"^ and recently formulated in 
a sophisticated fashion by Gull et al^^, which enabled a study 
at sufficiently low temperatures. 

In the QMC, we first expand the many-body partition func- 
tion in the interaction part of the Hamiltonian and apply the 
Hubbard-Stratonovich decoupling'*^ '*'* to the interaction part. 
This decomposes the many-body partition function into a sum 
of single-particle systems, which we collect. 

While the algorithm is based on the series expansion up 
to infinite order, it is feasible to obtain a numerically exact 
result because, after numerical convergence, all the orders in 
the expansion are virtually taken into account.'*" 

We implement the CTQMC sampling for the auxiliary 
fields, employing the updating algorithm proposed in Ref. 42. 
We typically take 2 x 10^ steps for each QMC simulation, 
and after convergence in the self-consistency loop, we aver- 
age over 30 data starting from go's which are different within 
a statistical error bar. 



III. RESULT AND DISCUSSION 
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FIG. 2. (color), (a) Low-energy structure of poles and zeros of G 
for t' = 0, U — 8t, and n — 0.91. Reproduction of Fig. 2(a) in 
Ref. 15. (b) A schematic picture for the relation between the cal- 
culated pseudogap and that observed in ARPES. The left side repre- 
sents the antinodal region around (tt, 0) while the right side the nodal 
region around (c) ReG~*^(k,cj) plotted against energy and 

momentum along the symmetric lines of k in the Brillouin zone. The 
parameters are the same as in (a). White curves between blue and red 
region represent poles of G while boundaries between black and yel- 
low region represent zeros, kpi and kF2 represent the Fermi wave 
numbers at the hole pocket along the nodal line, (d) Energy depen- 
dence of the self-energy at the Fermi momentum kF2 denoted in (c). 



A. Non-d-wave pseudogap 

First of all, we make a remark on the structure of pseudo- 
gap obtained in Fig. 2(a) in Ref. 15 [reproduced in Fig. 2(a)], 
which depicts the low-energy pole and zero surfaces calcu- 
lated by CDMFTh-ED for 9% hole doping to the Mott insu- 



lator at J7 = 8t and t' = 0. The pseudogap is formed by 
the zero surface (red) connecting two separated pole surfaces 
(green). Here we define the k-dependent amplitude A(k) 
of the pseudogap as the energy difference between the up- 
per and the lower poles at each k. Then we see in Figs. 2(a) 
and (c) that the gap in the direct transition opens in the en- 
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tire Brillouin zone, though it is somewhat larger in the antin- 
odal region [A((7r, 0)) = 0.35i] than in the nodal region 
[A(kFi = (0.557r,0.557r)) = 0.31i]. This is one of the cen- 
tral results of this paper, and clearly different from previous 
scenarios of pseudogap which assume a rf-wave gap above 
Tc, such as preformed pair,"* t/-density-wave,^''" resonating 
valence-bond, "'^^^ and nodal-liquid theory,'*^ '** since the di- 
rect transition gap closes in the nodal direction in these sce- 
narios. 



Nevertheless the electronic structure in Fig. 2(a), which was 
calculated without any assumption on the gap structure for the 
microscopic model, is consistent with the ARPES data, as de- 
scribed below. To start with, it should be noted that ARPES 
observes only the spectra below the Fermi energy i?F(= 0) 
if T is low, so that the gap amplitude is often estimated by 
symmetrizing the spectra below and above Eff^ This is a 
misleading, artificial procedure because it assumes a symmet- 
ric structure of the gap as to the Fermi energy and neglects 
the structure above i^p if any. Suppose we have the result in 
Fig. 2(a) only below iJp and symmetrize it as in the ARPES 
procedure to estimate the pseudogap, we end up with a rf-wave 
like gap since the gap in the part below is larger in the 
antinodal region while it is smaller or even zero in the nodal 
region [this is more clearly seen in Fig. 2(c) where we plot 
ReG^^ (k, a;) against energy and momentum along symmet- 
ric lines]. The situation is schematically shown in Fig. 2(b). 
Our results suggest that the "d-wave structure" is an artifact of 
the symmetrizing analysis, and in reality, the pseudogap has 
a non-rf-wave (namely, full-gap) structure. In APPENDIX C 
we confirm that the gap in the nodal direction persists for a 
larger cluster, by implementing an iVc = 8 CDMFTh-CTQMC 
calculation at a low temperature. 



It is interesting to examine this interpretation experimen- 
tally. Our theory predicts that a gap exists above E^ even in 
the nodal direction if the paramagnetic metal persists down 
to zero temperature. This means that the symmetrization"*^ 
used in ARPES breaks down in the pseudogap state due to the 
spectral asymmetry as to the Fermi energy. This may be ex- 
amined by spectroscopic or scattering probes to observe unoc- 
cupied electronic states above the Fermi level, such as the in- 
verse photoemission spectroscopy (IPES). It is highly desired 
to reveal the pseudogap structure by combining both PES and 
IPES without any symmetrization procedure. Other possible 
experimental probes suited for this purpose may be electron 
energy loss spectroscopy, resonant inelastic X-ray scattering, 
and time-resolved photoemission. In addition, another pos- 
sible study is ARPES on electron-doped cuprates. When we 
interpret the result with the electron-hole transformation, it 
provides information on the spectra of unoccupied states in a 
hole-doped system. We will discuss this in Sec. Ill I and III J, 
comparing the results with existing ARPES data on electron- 
doped cuprates. 




FIG. 3. (color online), (a) DOS for U = 8t, t' = 0, and 
n = 0.94 (solid curve) and for U = 8t, t' = -QAt, and n = 
0.95 (dashed curve), (b) Spatially averaged tunnelling conductance 
spectrum of lightly hole-doped cuprates, Nao.i2Cai.88Cu02Cl2 
(left scale) and Bi2Sr2CaCu208+a (right scale). Reproduction of 
Fig. 1(c) in Ref. 48. (c) Spectral function calculated with the 
CDMFT+ED around the antinode [(tt, 0)-(7r, f )] and the node 
[(^, ^)-(^, ^)] for t' = -QAt, U = and n = 0.95. The 
numerical data are broadened with -q = O.OSt. 



B. Spectral symmetry and asymmetry around the Fermi 
energy 

Although a thorough comparison of the calculated pole- 
zero structure with experiments is not possible at present be- 
cause of the lack of the experimental spectra above iJp, it is 
still significant to make comparison with available experimen- 
tal data which partly elucidated them. The scanning tunnelling 
microscopy (STM)"*** and the ARPES"^^ reported such data, 
where they found an electron-hole asymmetry of the spectra 
around E-p. The asymmetry provides a clue to the mechanism 
of the pseudogap, especially to the relation between the pre- 
formed pairing and the pseudogap, because the pairing will 
lead to a symmetric spectrum as estabUshed in the Bardeen- 
Cooper-Schrieffer (BCS) theory.^" 

Figure 3(a) shows the low-energy part of the density of 
states (DOS) for U = 8i, t' = 0, and n = 0.94 and for 
U = 8t, t' = -0.4i, and n = 0.95. The data show that 
the weight of the low-energy occupied states is significantly 
larger than that of unoccupied ones, in accord with other 
CDMFT studies. '■'•^^ The asymmetry is consistent with the 
spatially averaged electron-tunnelling spectra [Fig. 3(b), re- 
produced from Fig. 1(c) in Ref. 48] measured by the STM 
for lightly hole-doped cuprates, Nao.i2Cai 88CUO2CI2 and 
Bi2Sr2CaCu208+5, where the observed probability of elec- 
tron extraction is considerably greater than that of injec- 
tion. Moreover, the numerical data show a peak around 
u = 0.3t(0.4i) for t' = 0(-0.4t), which roughly agrees with 
the peak injection energies (100-300meV) seen in Fig. 3(b). 
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We note that the numerical data have the minimum above 
the Fermi level while the STM spectra have a V-shape gap 
with the minimum at zero bias. The V-shape gap might be 
attributed to the soft Coulomb gap*"' or soft Hubbard gap^^ 
caused by an interplay of electron correlations and random- 
ness, as indicated by the strong charge inhomogeneity ob- 
served by the STM at surfaces.^** 

ARPES allows a more detailed comparison of the spec- 
tra. Yang et al.*'^ succeeded in deriving low-energy (a; < 
O.OSeV) spectra above the Fermi level by carefully analyzing 
the ARPES data for Bi2Sr2CaCu208+5. This was done by us- 
ing the fact that ARPES data at finite temperatures contains 
information of unoccupied states because of the smeared tail 
of the Fermi distribution function for > 0. In the supercon- 
ducting state they observed electron-hole symmetric spectra 
both in the nodal and antinodal region, in accordance with the 
BCS spectral function. Meanwhile in the pseudogap state 
they found that (i) around the antinode the spectrum is nearly 
symmetric with intense peaks below and above Ep separated 
by a gap of ^ 0.06eV, and that (ii) as approaching the node the 
peak below Ep goes up and eventually crosses Ep while the 
peak above Ep disappears from the measured energy range. 
While (ii) clearly shows the asymmetry as to Ep, Yang et al.'^'^ 
interpreted (i) as an evidence of preformed pairing. 

Liebsch and Tong^^ obtained an asymmetry similar to (ii) 
with the CDMFT. Nevertheless it is still worthwhile to see 
if the symmetry (i) can be reproduced since it is relevant 
to the mechanism of the pseudogap. Figure 3(c) shows the 
normal-state spectral function yl(k, w), calculated with the 
CDMFTh-ED using 77 = OMt. We find that around the 
antinode, the two intense peaks reside nearly symmetric with 
opening of a gap. As approaching to nodal region, the gap 
below the Fermi level monotonically decreases and vanishes 
when the lower peak reaches the Fermi level [right panel of 
Fig. 3(c); see also Fig. 2(a)]. These behaviors are consistent 
with the experimental observations (i) and (ii). We note that 
in the t' — case [Fig. 2(c)] the gap at (tt, 0) above Ep is 
larger than that below Ep, whereas the spectrum is more sym- 
metric for t' = —OAt [Fig. 3(c)], which is more appropriate 
for BizSraCaCuaOs+a. 

In the nodal direction, the gap lies above w ~ O.lt. This 
is again consistent with the ARPES"*' which saw the region 
below O.OSeV ~ O.OSt and observed no gap in this direction. 
We note that this is distinct from the picture in Ref. 33, where 
the gap was assumed to close in the nodal direction. 

Thus we have shown that the ARPES data'*'' does not nec- 
essarily indicate a preformed pairing, but is rather naturally 
interpreted as a consequence of the dispersive zero surface. 



C. Back-bending behavior of dispersion 

Another important aspect captured in Fig. 2(a) is the back- 
bending behavior of the band cutting the Fermi level. Namely, 
the pole surface, which monotonically goes up to (vr, tt) in the 
bare dispersion, is bent back below Ep around (tt, tt). This 
is more clearly seen in Fig. 4(a), which is an enlarged view 
of Fig. 2(a) in a lower energy range. The back bending can 
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FIG. 4. (color online), (a) Enlarged view of Fig. 2(a) at low energy. 
Solid (orange) curves are guides to the eye for band dispersions in 
nodal and antinodal regions, (b) Binding energy of the ARPES spec- 
tral peaks along momentum cut with regular ky intervals. Repro- 
duced from Fig. 2(d) in Ref. 53. The data was taken for the pseudo- 
gap state (T = 140K) of a thin film Bi2Sr2CaCu208 sample with 
Tc = 90K. Ao(~ 50meV) is the maximum of the d-wave gap fitting 
the data. 



be seen below Ep around the antinodes while it is above Ep 
around the nodal direction. 

Actually the ARPES^-^ observed a similar behavior. Fig- 
ure 4(b) is a reproduction of Fig. 2(d) in Ref. 53, where the 
back-bending behavior is seen only around the antinode. In 
the light of Fig. 4(a) the absence of the back bending around 
the node will be simply because the top of the band is located 
at a higher energy than the maximum energy of the measure- 
ment. Note that i' < lifts (lowers) the band around the 
nodes (antinodes). Enhancing —t' shifts the Fermi surface to 
(0, 0) as well, as we will discuss in Sec. IIIH. Then a better 
agreement with the ARPES data may be reached. 

In Ref. 53 the back-bending behavior around the antinodes 
was interpreted as an evidence of preformed Cooper pairs in 
the pseudogap state because it resembles a band dispersion in 
the BCS superconductors.^*^ However, our result implies an- 
other simple interpretation: The band is pushed down by the 
neighboring zero surface, which cuts the Fermi level around 
(tt, tt), due to the large self-energy around it. This picture is 
confirmed in Fig. 2(d), where we plot the energy dependence 
of the self-energy at the Fermi momentum kpa closest to the 
zero surface. We see that the real part of the self-energy is 
negatively large below the zero at w ~ 0.06t, which should 
push down the band around the momentum. The result in- 
dicates that the pair formation is not necessary and the zeros 
of G resulting directly from strong correlation effects give an 
alternative interpretation of the back-bending dispersions. 

Recently, an evidence against the scenario that the back 
bending is a consequence of preformed pair was further 
reported. In this ARPES observation, the back-bending mo- 
mentum, i.e., the location of the pole at the lowest binding en- 
ergy along momentum cuts, is clearly deviated from the Fermi 
momentum at above the pseudogap opening temperature. In 
contrast, we note that the two momenta should agree in the 
BCS theory and in the preformed pair scenario as well. Fur- 
thermore, the ARPES reported that the back-bending momen- 
tum along (tt, 0)-(7r, tt) shifts closer to (tt, tt) than the Fermi 
momentum. This indicates a zero surface around (tt, tt) at the 
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FIG. 5. (a) Energy-distribution curves calculated with the 
CDMFT+ED along symmetric lines (0, 0)-(7r, 0)-(7r, tt) for t' = 
-~0.2t, U = 12t, n = 0.93, and 77 = Q.lt. For clarity the curves 
are offset by 0.06 (0.1) for (0, 0)-(7r, 0) [(tt, 0)-(7r, tt)]. (b) ARPES 
spectra for La2-iSri,Cu04 {x = 0.07), reproduced from Fig. 5 in 
Ref 56. 



Fermi level and a center of the gap residing above E^, in full 
consistency with our pole-zero structure Fig. 2(a). 

We note that a back-bending dispersion was already ob- 
served in an early QMC study,^^ where antiferromagnetic 
fluctuations were proposed as the origin of the pseudogap. Al- 
though our numerical data do not exclude the antiferromag- 
netic fluctuations from the possible mechanisms, the less-k- 
dependent pseudogap as well as the asymmetric location of 
the hole pocket (see Sec. IIIG) seems to oppose the mecha- 
nism; instead it rather supports that the pseudogap is a direct 
consequence of the proximity to the Mott insulator. 



(a) r =0, (;=i2f (b) r =-o.2t, u=^2t (c) r =-o.2t, u=m 

n=0.91 n=0.92 n=0.85 




(0,7l/8) (77t/8,7l) 



(d) r =0 ^=8? (e) ARPES for Bi2Sr2CaCu208+d 
n=0.91 [Graf et al., PRL 98, 067004 (2007)] 




D. Energy-distribution curve 

Zeros of G are not directly seen in spectra. Their footprints 
may, however, be detected through a sudden suppression of 
spectra due to a large ImS around the zeros. Figure 5(a) 
shows energy-distribution curves (EDC) of the spectral func- 
tion along momentum cuts (0, 0)-(7r, 0)-(7r, vr), calculated by 
the CDMFTh-ED with rj = O.lt for t' = -0.2t, U = 12t, and 
n — 0.93. The results exhibit a coherent peak around (tt, 0) 
just below the Fermi level, its shift to lower energy from (tt, 0) 
to (f-, 0), and the incoherent feature around (0, 0) and (tt, tt). 
All these features are consistent with ARPES data [Fig. 5 in 
Ref. 56, reproduced in Fig. 5(b)]. This agreement supports 
our zero mechanism and shows that the incoherent feature 
can be interpreted as the effect of zeros of G: Since the zero 
surface exists just above the band around (tt, tt) [Fig. 2(a)] 
and since another one is located just below the band around 
(0, 0) [not shown in Fig. 2(a), but can refer to Fig. 10(b) be- 
low], the spectrum associated with the band is smeared around 
these momenta due to the large ImS around the zeros. While 
the suppression around (tt, tt) is relevant to the emergence of 
Fermi arc,'^"'^'^" that around (0,0) is related to the waterfall 
behavior discussed in the next subsection. 



FIG. 6. (color online), (a)-(d) Spectral distributions calculated with 
the CDMFT+ED along the momentum cut (0, f )-(^,7r) [denoted 
by a red arrow in the inset of (a)] for a variety of parameter sets. The 
data are broadened with rj = 0.2t. Grey arrows denote the energies, 
El and E2, discussed in the text, (e) ARPES spectrum for overdoped 
Bi2Sr2CaCu208+i, reproduced from Fig. 1(c) in Ref. 58. 



E. Waterfall 

Recent ARPES^^"*' observed an anomalous spectral struc- 
ture, called "waterfall", at high binding energies ^ 0.3-0.4eV, 
where the band cutting the Fermi level suddenly loses the 
spectral weight and starts falling down to — 0.7eV with 
a suppressed intensity. Below the energy a strong intensity 
emerges again from nearly the same momentum as that the 
waterfall starts. This high-energy anomaly has been found 
quite generally in hole-doped cuprates, irrespective of the 
presence or absence of superconductivity, under, optimally or 
overdoped, and detailed compositions. Moreover similar fea- 
tures have been reported in other transition-metal compounds 
such as SVOs^^ and LaNiOa.*''' These suggest a universality 
of the phenomenon in strongly-correlated metals.*'* 

The dynamical cluster approximation (DCA) + QMC 
study^^ for the 2D Hubbard model with t' = reported a 
similar structure in the spectra. Based on the similarity of 
the DCA results to those with a perturbative calculation incor- 
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porating antiferromagnetic spin fluctuations, the authors pro- 
posed that the waterfall results from high-energy spin fluctu- 
ations. Here we implement a CDMFTh-ED study, and closely 
examine how the pole-zero structure underlies the waterfall 
phenomenon and how the energy scale of the waterfall de- 
pends on model parameters. Our analysis indicates alternative 
interpretation for the phenomenon. 

Figures 6(a)-(d) plot the spectral intensity calculated along 
the momentum cut from (0, ^) to tt) for various parame- 
ter sets. The results show nice resemblances with the ARPES 
result [Fig. 1(c) in Ref. 58, reproduced in Fig. 6(e)], i.e., (i) 
an abrupt change in the slope of the band around u = Ei, 
accompanied by the simultaneous reduction of the intensity, 
(ii) nearly vertical dispersion with suppressed intensity in the 
waterfall region, and (iii) reemergence of a strong intensity 
around uj — E2- The energy scale also agrees roughly with 
the experimental value, if one keeps in mind t 0.4 eV. 

We investigate how Ei and E2 depend on t', n, and U. 
First, comparing Figs. 6(a) and (b), we see that t' < reduces 
both \Ei \ and \E2\. Second, comparison of Figs. 6(b) and (c) 
shows that doping also reduces the energies. This doping de- 
pendence is qualitatively consistent with the DCA results^^ for 
16 sites cluster at a finite temperature. Third, the decrease of 
U from 12t [Fig. 6(a)] to 8t [Fig. 6(d)] increases the energies. 

These parameter dependences can be naturally explained as 
follows. In general, strong correlation effects make the width 
of coherent band (i.e., the band above Ei) narrower by in- 
creasing the mass of the low-energy particles. This leads to 
the decrease of |£'i| as U increases. Indeed the slope of the 
coherent band decreases from 1.3ta (a: the lattice constant) 
to I.Ota as U increases from 8t to 12t. Meanwhile t' lifts 
the band around the nodal direction, as will be discussed in 
Sec. IIIH. This elevates Ei and E2- The doping dependence 
seen in Figs. 6(b) and (c) can be qualitatively understood as 
a downward shift of the chemical potential with doping. As 
we showed in Ref. 15, in hole-doped Mott insulators doping 
causes a rigid-band-like shift of the chemical potential, though 
it cannot be described within the single-electron picture. In 
fact, the low-doping phase is not adiabatically connected with 
the Fermi liquid phase at higher dopings because of the inter- 
vening Lifshitz transition and the zero-surface emergence. 

The above interpretation suggests that the energy scale of 
the waterfall is not necessarily related to the antiferromagnetic 
spin fluctuations^^ but is considered to be a rather direct con- 
sequence of electron correlation, i.e., the mass renormaliza- 
tion of the coherent band. The picture is corroborated by the 
fact that a similar behavior can be seen within the single-site 
DMFT.''4 

In view of the pole-zero structure in Fig. 1(b) in Ref. 15, 
the waterfall emerges in the energy region —2t ^ w < —t, 
where many smeared pole and zero surfaces pile. The con- 
gestion of poles and zeros within the finite-cluster calculation 
indicates an incoherent nature of this energy-momentum re- 
gion, which would result in the waterfall. Notice that the wa- 
terfall obtained in this paper has neither a gap nor dispersive 
features^**, but a vertical structure with a suppressed inten- 
sity. This unusual structure comes out thanks to the use of the 
cumulant periodization (see Sec. II) which can describe the 
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FIG. 7. (color online). Spectral function calculated with the 
CDMFT+ED for t' = 0, U ^ 8t, n = 0.93, and 77 = O.OOSt 
along the momentum cuts (a) (0, 0)-(7r, vr) and (b) (tt, 0)-(7r, tt). The 
(blue) dashed lines in (a) are guides to the eye. 



strong momentum dependence of the self-energy. Below the 
pile we see a relatively coherent pole surface extending down 
to ^ — 3.5i. This corresponds to the band reemerging at the 
high binding energy [(iii)]. Note that a comparison at higher 
binding energies is difficult because the ARPES spectra are 
overlapped by other Cu-d or O-p bands. 



F. Low-energy kink in dispersion 

In prior to the discovery of the waterfall, another anomaly 
(kink) in dispersion has been found in ARPES at lower bind- 
ing energies^^ and its origin has been in dispute. The kink 
is located at a binding energy around O.OSeV, where the band 
cutting the Fermi level sharply changes its slope. In the nor- 
mal state the kink is clearly seen around the nodal region 
while, as approaching the antinodal region, the band becomes 
flatter and then the kink becomes less visible. The band struc- 
ture continuously changes from the nodal to antinodal region 
where a pseudogap exists just above the band.^^'^^ 

The CDMFTh-ED results in Fig. 7 show similar behaviors: 
(i) The band suddenly changes the slope at a kink around uj — 
— O.lt ^ — 0.04eV in the nodal region [Fig. 7(a)], and (ii) as 
approaching the antinode, the band becomes flatter and the 
slope change becomes weaker [Fig. 7(b)]. 

In general, the zero surface which generates the pseudogap 
pushes down the dispersion near the Fermi level. This makes 
a quick change of the slope of the dispersion distinct from the 
part deeply below the Fermi level, which may be the underly- 
ing origin of the kink formation. However, the sudden change 
of the dispersion observed in Fig. 7(a) naturally requires a pre- 
cursory formation of a tiny zero surface around uj — —O.lt. 
This means a coupling of the quasiparticle with some other 
excitations with this energy, whichever bosonic or fermionic. 
The energy resolution of the present cluster size is not obvious 
and it could be an artifact of the present small cluster calcula- 
tion. Nevertheless, it is remarkable to note that the observed 
kink in Fig. 7(a) is rather universal and we see similar kinks 
in other part of the quasiparticle dispersions, for instance in 
Fig. 10(b) below. This may be alternatively interpreted that. 
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FIG. 8. (color online). Momentum maps of (a) ReG ^(k, = 0) 
and (b) A(k,uj = 0) calculated with the CDMFT+ED for t' = 0, 
U = 12t, and n = 0.91. 77 = O.OOOU (O.Olt) is used for (a) [(b)]. 
CDMFT+CTQMC results of -/3G(k, |) for t' = 0, U = I2t and 
T = 0.025t at (c) n = 0.95 and (d) n = 0.90. (e) Momentum- 
distribution curves calculated with the CDMFT+ED along (0, 0)- 
(tt, tt) for the same parameters as those for (a) and (b). rj = O.Olt 
(solid curve) and O.lt (dashed curve). 



in the strongly correlated region, the quasiparticle is strongly 
renormalized and may couple to various intrinsic electronic 
modes of charge and spin origins. 

In high-Tc cuprates the mechanism of the kink has been 
extensively discussed in the literature from the viewpoint of 
a coupling of electrons with some bosonic mode, such as 
phonon'''^ and magnetic ones.*^'*'^ Our result, however, implies 
that the kink may be a rather general phenomenon in the prox- 
imity to the Mott insulator,^"* while the specific mechanism of 
kink formation is left for future studies. 



G. Fermi arc and hole-pocket Fermi surface 

Next we shift our focus on the zero-energy electronic struc- 
ture, namely, Fermi surface. The spectra in hole-doped 
cuprates have been extensively studied by ARPES.^ One of 
the most remarkable findings in these studies is the observa- 
tion of truncated Fermi surfaces, called "Fermi arc".'*^ 

Several authors have discussed the Fermi arc in terms of 
coexisting poles and zeros of G. i3, 15,30,32-36 jjj jj^jg mechanism 
hole-pocket Fermi surfaces around the nodal directions coex- 
ist with a zero surface around (7r,7r). Because ImS is large 
around the zero surface, the pockets lose the spectral intensity 
more on the side closer to the zero surface, leaving arc-like 



spectra on the opposite side. This is reproduced in Figs. 8(a) 
and (b). Figure 8(a) shows the underlying pole-zero structure 
calculated by the CDMFTh-ED at T = without any substan- 
tial smearing to the singularities, where a clear hole-pocket 
Fermi surface coexists with a zero surface. When we intro- 
duce a smearing to the singularities (see Sec. II A), we obtain 
the spectral map. Fig. 8(b), in which the arc-like structure is 
observed. 

We confirm the zero mechanism by directly calculating 
a spectral weight at finite temperatures by employing the 
CTQMC method as the solver for the CDMFT. To obtain the 
spectrum without through any analytic continuation procedure 
which inevitably suffers from a large error bar, we calculate 



-/3G(k,r = /3/2) = | r 



A{-k,i 



cosh(/3w/2) 



dio 



(9) 



as a function of k. This quantity is an integral of the spec- 
tral weight over a width ~ T around u = and approaches 
yl(k, ll! — 0) for large (3, so that gives an estimate of A(k, uj — 
0) at low temperatures.^" In Figs. 8(c) and (d) we plot the 
resuks at T = 0.025t. For n = 0.95 [Fig. 8(c)] we see 
an arc at a location similar to the one in Fig. 8(b) while for 
n — 0.9 [Fig. 8(d)] the spectra extend to the antinodal regions. 
The doping evolution of the spectra is qualitatively consistent 
with that obtained with the CDMFTh-ED in Fig. 2(c)-(e) in 
Ref. 15. The qualitative agreement between the CDMFTh-ED 
result with broadening 77 and the CDMFT+CTQMC resuk at 
finite temperatures supports that the phenomenological smear- 
ing factor rj well simulates the thermal effects. Moreover 
it corroborates the above-mentioned zero mechanism for the 
emergence of the Fermi arc. We note that the zero surface 
around (tt, tt) is also consistent with the results by other clus- 
ter schemes such as DCA, which observed a strong scatter- 
ing amplitude in the momentum patch around (tt. tt) in both 
A^c = 4 (Ref. 71) and A^^ = 8 (Ref. 29) calculations. 

Interestingly, recent high-resolution ARPES^^ reported the 
existence of hole-pocket Fermi surfaces around the nodal di- 
rections in underdoped Bi2Sr2-a;Laa;Cu06+5. The observed 
hole pockets have much less intensity on the side closer to 
(tt, tt) than the opposite. This is consistent with the zero 
mechanism in Fig. 8(b). Important findings in the ARPES 
are that (i) the hole pockets are not located symmetrically 
with respect to the antiferromagnetic Brillouin zone boundary, 
i.e., (tt, 0)-(0, tt) line (and its symmetrically-related ones), 
and that (ii) the spectral intensity on the (tt, tt) side is finite 
even in the nodal directions, (i) excludes several scenarios 
for the pseudogap, which attribute the pockets to symmetry 
breakings^'^ '" because in these theories hole pockets should 
be centered symmetrically as to the (tt, 0)-(0, tt) line.^^ Mean- 
while (i) is consistent with our zero mechanism that does not 
assume any symmetry breaking, (ii) is at odds with the theory 
based on an ansatz that gives a zero intensity in the outer point 
of the pocket crossing the nodal direction,^^ while it agrees 
well with our numerical data in Fig. 8(b). This is more clearly 
seen in Fig. 8(e), where we plot A(k, 0) along the momentum 
cut in the nodal direction: For ry = O.Olt, in addition to the 
main peak at k ~ (O.SBtt, O.SGtt), we see the secondary peak 
at ~ (O.Ttt, O.Ttt), corresponding to the hole pocket structure. 
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FIG. 9. (color online), t' dependence of -/3G(k, f ) in Eq. (9) at 
T = 0.05t for (a)(b)(c) U ^ Stand n = 0.93, and (d)(e)(f) U = 12t 
and n — 0.90. White dashed line represents the antiferromagnetic 
Brillouin zone boundary. 



The second peak is, however, not visible for r] = O.lt and 
only the broadened main peak can be seen there. Namely, 
in our theory, the spectrum looks either a pocket like or an arc 
like depending on how the incoherence due to the zero surface 
is strong. The strength of the incoherence is controlled by the 
value of J] or temperature, or the distance between zeros and 
poles, as partly demonstrated in Fig. 4 in Ref. 15. Figure 8(e) 
also suggests that the energy resolution required to detect the 
pocket in the example of the pole-zero structure as Fig. 8(a) 
is in between 0.01t{^ 4meV) and 0.1i(~ 40meV), which 
roughly corresponds to the energy resolution in ARPES. This 
consideration on the incoherence explains why the pocket had 
not been detected until the recent high-resolution ARPES. 

One obvious difference between Fig. 8(b) and the ARPES 
is the location of the pocket: The pocket in Fig. 8(b) resides 
closer to (tt, tt) than that observed in the ARPES. This may 
be attributed to the difference between the model parameters 
we used and realistic ones for Bi2Sr2-2:La2:Cu06+5. In par- 
ticular, we show in the next subsection that the next-nearest- 
neighbor transfer t', which is zero in Fig. 8(b), indeed shifts 
the pocket in the direction to (0, 0). 



(olto 
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(c) t =0, n=^.^2 



(d)r =-0.2t, n=1.13 
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FIG. 10. (color). ReG~^(k,a;) calculated with the CDMFT-^ED 
along symmetric lines for U — 12t. (a) t' — and n = 0.91, 
(b) t' = -0.2t and n = 0.92, (c) = and n = 1.12, and (d) 
t' — ~0.2t and n — 1.13. 77 = O.OOlt is used. Color scale is the 
same as that in Fig. 2(c). 



alistic determination of the model parameters, longer-range 
transfer integrals and Coulomb interactions, and a calculation 
on a larger cluster. These remain for future researches. 

The above shift of the arc with t' can be understood by plot- 
ting the underlying pole-zero structures. Figures 10(a) and (b) 
compare the structures for t' = and t' = — 0.2< in a lightly 
hole-doped region. The figures show that at low energies t' 
lifts the poles in the nodal direction while lowers them around 
the antinodes, as already expected from the bare dispersion 
Eq. (2). Then the hole pocket expands and the Fermi surface 
with the stronger intensity shifts to (0, 0) direction. We note 
that, as we discuss in the next subsection, the third-neighbor 
transfer integral, which is not taken into account in the present 
calculation but exists in real materials, further enhances the 
pocket around 



H. Effect of the next nearest-neighbor transfer t' 

Here we systematically study the effect of t' on the elec- 
tronic structure in lightly hole- and electron-doped regions. 
Figure 9 shows t' dependence of the integrated spectra, 
Eq. (9), calculated with CDMFTh-CTQMC at T = O.G5t for 
(a)-(c) U ^ 8t and n ^ 0.93 and (d)-(f) U = I2t and 
n = 0.90. We see that the arc shifts to (0, 0) as -t' /t in- 
creases, and that at t' = —0.4^, which is a reasonable value 
for Bi2Sr2-j;La^Cu06+5, the arc resides inside of the antifer- 
romagnetic Brillouin zone, in consistency with experiments.^^ 
We also see that t' increases the curvature of the arc. A com- 
parison with ARPES in further detail will require a more re- 



The effect of t' is totally different for electron-doped cases. 
Figures 10(c) and (d) show the structures for t' = and —0.2t, 
respectively, in an electron-doped region. The effect of t' at 
low energies is again understood by Eq. (2): t' lifts the poles in 
the nodal direction while lowers around the antinodes. How- 
ever, the consequence on the Fermi surface completely differs 
from that in hole-doped cases. For the particle-hole symmet- 
ric t' = case, the electron pockets appear around the nodes 
for electron doping, corresponding to the hole pockets in hole- 
doped cases, however, for finite t', the electron pockets appear 
around the antinodes, as shown in Fig. 10(d). This is because 
the low-energy zero surface pushes up the dispersion around 
(0, -1). We compare the Fermi surface structure with ARPES 
for the electron-doped cuprates in detail in Sec. Ill J. 
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FIG. 11. (color online), (a) [(b)] DOS for t' = -0.2t, U = I2t 
and n = 0.92(1.09) with r? = O.Olt. LHB and UHB represent the 
lower and upper Hubbard bands, respectively. (c)(d) ARPES spectra 
and (e)(f) their second derivatives in uj for Ndi.85Ceo.i5Cu04. Re- 
produced from Fig. 1 in Ref. 73. (g)(h) Spectra calculated with the 
CDMFT+ED for t' = ~0.4t, U = 12t, and n = 1.11 along mo- 
mentum cuts, (0, 0)-(7r, tt) and (0, Tr)-{n, n), respectively. 77 = 0.3t 
is used. 



I. High-energy spectra in electron-doped cuprates 

Recently high-energy spectra of an electron-doped cuprate 
Ndi 85Ceo.i5Cu04 were studied by ARPES. To discuss the 
results we first emphasize that here the ARPES observes 
a totally different energy region from that for hole-doped 
cuprates. This is illustrated in Figs. 11(a) and (b): ARPES 
observes the lower Hubbard band (LHB) in hole-doped cases 
while in electron-doped cases it does ingap states between the 
Mott gap and the pseudogap, in addition to the pseudogap 
structure itself. As mentioned in Sec. IIIB, with the electron- 
hole transformation the ARPES result can be interpreted as 
that for unoccupied states in a hole-doped system. Therefore 
the results in Ref. 73 provide a precious opportunity to com- 
pare the calculated spectra of ingap states with experiments. 
The pseudogap and the Fermi surface are discussed in the next 
section. 

The ARPES spectra and its second derivative in uj (Fig. 1 
in Ref. 73) are reproduced in Figs. 1 l(c)(d) and (e)(f), respec- 
tively. Around the antinode the spectra [Fig. 11(d)] suddenly 
lose the intensity at ~ — 0.4eV. The analysis on the sec- 



ond derivative [Fig. 11(f)] shows that the low-energy band 
reaches the bottom at this energy. On the other hand, around 
the node the low-energy band persists down to w ~ — 0.7eV 
[Fig. ll(c)(e)]. 

Figures 1 1(g) and (h) depict the spectra calculated with the 
CDMFTh-ED for t' = -0.4t, U = Ut, and n = 1.11 along 
momentum cuts, (0, 0)-(7r, vr) and (0, tt)-{tt, tt), respectively. 
For a; < we see that the band bottom is located at a; ~ — 0.7t 
around (0,7r) while it is located at a higher binding energy 
a; ~ — i around (0,0). This is consistent with the ARPES 
data. Note that the numerical data do not show any inten- 
sity below the band bottoms, w < —0.7t {—t) around (0, tt) 
[(0, 0)], which is within the gap region between LHB and the 
ingap states, whereas the ARPES spectra show waterfall-like 
structures in — leV < w < — 0.7eV (lu < — 0.4eV) around 
(0, vr) [(0, 0)]. The origin of this discrepancy is not clear, but 
one possibility is a multiband effect, which is not considered 
in the present single-band model. 

Meanwhile for w > 0, we see a waterfall in the dispersion 
at w ~ t-2t. This corresponds, with the electron-hole trans- 
formation, to the waterfall in hole-doped cases. The result 
predicts that this anomaly will be observed in electron-doped 
cuprates when the unoccupied spectra become available up to 
the high energy.^'* 



J. Pseudogap and Fermi pockets in electron-doped cuprates 

Lastly we discuss low-energy electronic structures in 
electron-doped cases. The pseudogap has been experimen- 
tally observed also in electron-doped cuprates, although 
the doping range is much more limited than that in hole- 
doped cases due to the wider antiferromagnetic region near 
the undoped Mott insulator. The mechanism of the pseu- 
dogap has extensively been discussed within weak-coupling 
theories^^'^^, which has successfully reproduced various ex- 
perimental results such as the doping evolution of the Fermi 
surface. In these theories the mechanism has been ascribed 
to the antiferromagnetic long-range correlations. On the other 
hand, in the strong-coupling regime the cluster-perturbation 
theory^'^ and the CDMFT^''^^ found that the pseudogap re- 
sults from nonlocal but short-ranged dynamics without any 
long-range order or correlations. As pointed out in Ref. 20 
and 22, the gap amplitude does not scale as J ^ /U , so 
that the short-ranged dynamics is not simply ascribed to the 
AF correlation. Because it is unknown which mechanism, 
weak-coupling or strong-coupling, is relevant to electron- 
doped cuprates, it is worthwhile to see whether and how ex- 
perimental data are understood in the strong-couphng picture, 
especially from the perspective of zeros of G. 

The available ARPES data for the electron doped cuprates 
seems to have relatively poor resolutions and a clear dis- 
persion has not been reported. Nevertheless, a recent 
ARPES^^ measurement appears to have revealed the struc- 
ture of the pseudogap for a family of electron-doped cuprates, 
Ln2-xCtxCviOi {Ln = Nd, Sm, and Eu). In Fig. 2 in 
Ref. 78 the EDC peak position jumps from oj ~ — O.OSeV 
to u! ^ — 0.15eV around antinode, which implies the pres- 
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(0,0) (7t,0) 




FIG. 12. (color online), (a)-(c) ARPES spectra for Nd2_ a; Ce;j,Cu04. 
Reproductions of Fig. 3 in Ref. 79. (d)-(f) -^G(k, f ) calcu- 
lated witii tlie CDMFT+CTQMC for t' = ~QAt, U = 12t, and 
T = 0.033t at n = 1.05, 1.1, and 1.15, respectively, (g)-(i) 
ReG"^ (k, uj = 0) calculated with the CDMFT+ED for t' = ~0.4t, 
U =12tatn = 1.14, 1.17, and 1.19, respectively, r; = O.OOOlt is 
used. 



ence of the pseudogap opening below the electron pocket. 
This highly asymmetric position of the pseudogap is a feature 
which cannot be seen by ARPES for the hole-doped cuprates 
but observable for electron-doped ones: Because in electron- 
doped cases the zero surface forming the pseudogap mainly 
extends below Ep [see Fig. 10(d)], the ARPES can observe 
the major part of the pseudogap, in contrast to the hole-doped 
cases. Around the node, the EDC peak position saturates at 
UJ ~ — O.OSeV for Ln =Eu compound, implying the pseu- 
dogap at a; > —O.OSeV, while the pseudogap is not clear for 
Ln =Nd in EDC. This is consistent with the observation in 
Fig. 10 that the pseudogap position in the nodal region is low- 
ered for smaller \t' /t\ and that the Eu compound appears to 
have a relatively small \t' /t\. Namely, the main part of the 
pseudogap in the nodal region of the Nd compound is located 
in the positive energy side, while it becomes visible below 
Ep for the Eu compound. This suggests a gap around the 
node, which can be interpreted as an indirect evidence of a gap 
opening above Ep in hole-doped cases via the electron-hole 
transformation. These observations are qualitatively consis- 
tent with the pole-zero structure in Fig. 10(d), supporting the 
non-rf-wave (fully-opened) pseudogap proposed in Sec. Ill A. 
We propose to perform ARPES measurements of the electron 
doped compounds with a better resolution than Ref. 78, be- 
cause the pseudogap appears to be suggested only by the jump 
in EDC and the detailed pseudogap structure in the momen- 
tum space is not very clear so far. 

The ARPES^**-^"^ has also revealed a characteristic evolu- 



tion of the Fermi surface with doping. Figures 12(a)-(c) are 
the reproductions of the ARPES results (Fig. 3 in Ref. 79) on 
Nd2-,cCe:rCu04. The ARPES study found that (i) at small 
doping (x = 0.04) strong intensity emerges around (7r,0) 
and its symmetrically-related points [Fig. 12(a)], (ii) as dop- 
ing is increased, additional intensity develops around (§, §) 
[Fig. 12(b)], and (iii) at large doping the spectra merge, form- 
ing a single large surface around (tt, tt) [Fig. 12(c)]. In Ref. 79 
it is also seen that the band cutting the Fermi level at around 
evolves with doping from a high to low binding energy 
while that around (tt, 0) emerges first around the Fermi level 
and extends to a high binding energy. This suggests that the 
Fermi surface around ( f , ) is a hole pocket while the one 
around (tt, 0) is an electron pocket. 

The above features (i) and (iii) were qualitatively repro- 
duced by the cluster perturbation theory^" and by CDMFT 
with the E^' and G^^ periodizations on the 2D Hubbard model 
at T = 0. In Ref. 20 and 21 the suppression of the intensity 
around the node at low dopings was understood by the pres- 
ence of the hot spot, where ImE is large, in this region. 

In the following we clarify how the above observations (i)- 
(iii) can be understood in terms of the underlying zero surface. 
First, we show the CDMFTh-CTQMC results at T = O.Obt. 
Figures 12(d)-(f) depict the integrated low-energy spectra, 
Eq. (9), for t' = -0.4^ U = 12i, at 5, 10, and 15% electron 
dopings, respectively. At n = 1.05 a strong intensity emerges 
only around the antinodal points, in accord with (i). A large 
Fermi surface around (vr, tt) at n = 1.15 agrees with (iii). As 
to (ii), our results do not show strong intensity around 
this will be attributed to the absence of the third-neighbor 
transfer integral t" in the present calculation, as discussed be- 
low. 

Second, Figs. 12(g)-(i) show the CDMFTh-ED results for 
t' = -OAt and C/ = 12t at n = 1.14, 1.17, and 1.19, re- 
spectively. We see that a zero surface exists at the Fermi level, 
surrounding (0,0). For low dopings electron-pocket Fermi 
surfaces reside around the antinodal points. As the doping is 
increased, the pockets expand. Then they merge around the 
nodal point and change the topology into two Fermi surfaces, 
one around (7r,7r) and the other around (0,0). Further dop- 
ing makes the inner Fermi surface annihilate in pair with the 
zero surface, resulting in the Fermi liquid with the single large 
Fermi surface around (vr, tt).***' Thus in electron doping there 
occur, at least, two phase transitions, the Lifshitz transition**' 
and a pole-zero annihilation transition, on the way from the 
Mott insulator to the Fermi liquid. 

The presence of the zero surface accounts for the weak 
spectral intensity around the nodal direction in Figs. 12(d) and 
(e). Then how is the intensity around the node in Fig. 12(b) 
and (c) explained? To answer this we consider how the elec- 
tronic structure. Fig. 10(d), changes with t", which is not in- 
cluded in our calculations based on the 2x2 cluster Although 
\t"\ is usually several times smaller than \t'\ and therefore 
negligible for discussing global electronic structures, it can 
be crucial to understand the nodal intensity. This is because 
around (f, f ) the pole surface at ~ — O.lSt in Fig. 10(d) 
moves up with —t'/t, as seen in comparison with Fig. 10(c), 
and is located just below Ep for —t'/t — 0.4 (not shown). 
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Since to the first approximation t"{> 0) elevates dispersions 
around (|, ^) while lowers around (0, 0), (tt, tt), (tt, 0), and 
(0, tt), according to 

e(k) ^ e(k) - 2t"[cos(2/c^) + cos(2fcj^)], (10) 

a hole pocket can emerge around (f , f )• This hole pocket 
will show up just inside the zero surface in Figs. 12(g) and (h). 
Recent experimental observation^^ of the Shubnikov-de Haas 
oscillation in Nd2-a;Cea;Cu04 also indicates a presence of the 
hole pocket around optimal doping. Therefore, to confirm the 
above scenario with a larger cluster calculation incorporating 
t" is an intriguing future problem. 

We note some additional indications of the zero surface 
in experiments. In Fig. 2(c)-(e) in Ref. 79, which plotted 
the ARPES EDCs for Nd2_a;Ce2:Cu04 along the underlying 
Fermi surface, a clear jump of the peak positions can be seen 
between the nodal and antinodal directions. This indicates the 
presence of the zero surface around the jump, in accord with 
the above picture [see Fig. 10(d)]. In Fig. 1 in Ref. 78, which 
plotted the ARPES spectra at the Fermi level, a similar struc- 
ture to Fig. 12(c) was found while it clearly shows a finite in- 
tensity in some regions on the (tt, 0)-(0, tt) line. This appears 
to be inconsistent with the zero surface on the (tt, 0)-(0, tt) 
line assumed in Ref. 33 for hole-doped cases, but compatible 
with our zero surface around (0, 0) as seen in Figs. 12(g) and 
(h). 



IV. SUMMARY AND CONCLUSION 

In summary, we have shown that various spectral anoma- 
lies observed in the pseudogap states of hole- or electron- 
doped cuprates are naturally understood in terms of under- 
lying pole-zero structure of electronic Green's function. The 
pole-zero structure has been calculated for the paramagnetic 
metallic phase in the 2D Hubbard model with employing the 
CDMFT-i-ED at T = and the CDMFT-hCTQMC at T > 0. 
We have confirmed that the cumulant periodization scheme^^ 
is suited for the parameter region of our interest, namely 
lightly doped Mott insulator in the intermediate to strong cou- 
pling region, since the cumulant is well localized within the 
cluster that we employed, as demonstrated in APPENDIX B. 
This corroborates the convergence and reliabihty of our mo- 
mentum resolution in the studies on Fermi surface topology 
and the differentiation in the momentum space. Furthermore, 
the cumulant periodization is useful since it enables to study 
the coexistence of pole and zero surfaces. 

The result calculated by CDMFTh-ED shows a pseudo- 
gap in the lightly doped region around the Mott insulator. 
The pseudogap is characterized by a low-energy zero surface, 
which connects the bifurcated low-energy bands [Fig. 2(a)]. 
In hole-doped cases the lower band cuts the Fermi level while 
the upper one resides at the energy higher than and adjacent to 
the pseudogap above the Fermi level. This directly leads to (i) 
the non-rf-wave character of the fully opened pseudogap and 
(ii) the spectral asymmetry around the Fermi level. Looking 
into the structure around (tt, tt) we have also found that a large 



ReE around the zero surface causes (iii) the back-bending dis- 
persion. All of (i), (ii), and (iii) are consistent with exper- 
iments semiquantitatively. The agreement imposes a strong 
constraint on theories for pseudogap: Although (ii) and (iii) 
have been interpreted as evidences of the preformed pairing 
in the literatures, (i) is in sharp contrast to d-wave-gap scenar- 
ios including those by preformed pair. It is compatible neither 
with other precursory or real rf-wave-type symmetry break- 
ings such as d-density wave nor with commensurate antifer- 
romagnetism. We have confirmed a full gap formation at the 
nodal point in a larger-cluster calculation in APPENDIX C. 
Moreover, the non-d-wave gap is supported by the agreement 
of our result on the pseudogap structure with that observed by 
ARPES for electron-doped cuprates. 

The zero mechanism also enables a simple and unified 
understanding of various spectral anomalies: In hole-doped 
cases the same zero surface which causes the back-bending 
behavior around the antinode induces the incoherence around 
(tt, tt), Fermi pockets and Fermi arcs, while a pile of pole and 
zero surfaces at a higher binding energy results in the high- 
energy kink (waterfall) and the incoherence around (0,0). 
Moreover we have found a low-energy kink structure in the 
nodal dispersion. On the other hand, in electron-doped cases 
the zero surface which constitutes the pseudogap crosses the 
Fermi level around (0,0), making electron pockets around the 
antinodes. All these features are consistent with experimen- 
tal results. We would like to emphasize that the zero surface 
does not result from symmetry breakings, but is a direct con- 
sequence of the strong correlation effect, i.e., proximity to the 
Mott insulator. 

To elucidate the effect of zeros at finite temperatures, we 
have also implemented CDMFTh-CTQMC calculations. The 
results qualitatively agree with those of CDMFTh-ED with a 
smearing factor 77, which confirms that the thermal scatterings 
broaden the effect of zeros and indeed create the characteristic 
spectra similar to those observed in the cuprates. 

To conclude, the origins of various anomalies in the elec- 
tronic structure of the normal state in the high-Tc cuprates 
are unified into the presence of the low-energy zero surface 
which persists against doping. The zero surface is also ex- 
pected to cause anomalous metalUc behaviors in other physi- 
cal quantities, such as the specific heat, electronic resistivity, 
and Hall coefficient.' Comparisons of these quantities with 
experiments remain for future studies. 

Although the microscopic mechanism creating the zero sur- 
face has not been discussed in this paper, it has emerged 
clearly as the proximity to the Mott insulator. The similarity in 
the structure of the zero surfaces between the pseudogap and 
the Mott gap in the undoped system [see Fig. 1(a) in Ref. 15] 
also implies the Mott origin of the pseudogap. At the same 
time, however, the pseudogap does not appear as that directly 
or continuously connected to the zero surface which forms 
the Mott gap in the undoped state. This is because, for exam- 
ple in the hole-doped case, the zero surface creating the fully- 
opened pseudogap is bounded above by the existence of low- 
energy excitations far below the upper Hubbard band (UHB). 
The reason why the ingap states are separated by the pseudo- 
gap from the main quasiparticle band is left for an important 
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future issue. In Ref. 15 we proposed a scenario in which the 
emergence of such a zero surface in the low-energy scale and 
the resultant pseudogap are a remnant of the binding of dou- 
bly occupied site (doublon) and empty site (holon) drastically 
weakened by the screening of Coulomb repulsions by mobile 
carriers. This has something to do neither with any symmetry 
breaking nor with its precursor Recently, it has been proposed 
that a gap arising from hybridization between the quasiparticle 
and a composite fermion excitation is responsible for the pseu- 
dogap and the related zero surface. Possible supplementary 
role of strong antiferromagnetic fluctuations are left for future 
studies. Further study on the mechanism of pseudogap will be 
described elsewhere. 
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APPENDIX A: BREAKDOWN OF S PERIODIZATION FOR 
THE MOTT INSULATOR 

Here we note that a standard periodization technique, which 
substitutes Q = S to Eq. (3), cannot reproduce the Mott 
gap for large U . Because ImG = in the Mott gap, ImE 
must be or oo in the whole Brillouin zone. However, this 
situation does never occur in Eq. (3) with Q = S unless 
e CjImE^j = or = oo for ah uj inside 

the gap, and these conditions are never satisfied as far as both 
t and U are nonzero and finite. 

As a matter of fact, we see the electronic structure shown 
in Fig. 13(a) when we use the S-periodization procedure for 
the same Mott insulator shown in Fig. 1(a) in Ref. 15 [and 
reproduced in Fig. 13(b) for comparison]. We see that zeros 
of G exist only at the Fermi level without a dispersion, and 
that poles of G extend to the Fermi level around (0, 0) from 
the positive cj side and around (tt, tt) from the negative cj side, 
making the density of states half metallic. This failure of the 
E periodization in the Mott insulator is ascribed to the fact that 
E is not localized within the 2x2 cluster The nonlocality of 
S is a direct consequence of the presence of dispersive zeros 
of G, i.e., momentum-dependent divergence of E. 

As we have discussed so far, zeros of G still persist in doped 
Mott insulators up to a critical doping level beyond which the 
Fermi liquid emerges. Therefore S should be highly nonlocal 
also in the non-Fermi-liquid region and the S periodization 
again breaks down there. On the contrary, the cumulant M 
is well localized within the 2x2 cluster in this region, as 
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FIG. 13. (color online). Re*:?"^ (k, lj) obtained with (a) E and (b) 
M periodization for a parameter region of the Mott insulator {t' = 0, 
U = 8t, and n = 1). (c)(d) Density of states corresponding to (a) 
and (b), respectively. 
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FIG. 14. (color online). (a)(b) Matsubara-frequency dependence of 
the real and imaginary parts of cumulants M for various real-space 
vectors r = (i, j) in the 4x4 cluster [shown in the inset of (b)], 
calculated with the CDMFT+QMC for t' ^ Q, U ^ St, n ^ 0.91, 
and T = Q.lt. (c) The cumulants at the lowest Matsubara frequency 
plotted against the Euclidean distance. 

we mentioned in Sec. II. In APPENDIX B we give another 
numerical evidence for the locality of the cumulant in doped 
Mott insulators. 



APPENDIX B: LOCALITY OF CUMULANT IN DOPED 
MOTT INSULATORS 

Here we present a CDMFTh-QMC result for an iVc 4 x 4 
cluster in a parameter region of doped Mott insulators. We use 
the Hirsch-Fye algorithm^^ and calculate the cluster cumulant 
M'^ for i' = 0, [/ = 8t, n = 0.91, and T = Q.lt. 

Figures 14(a) and (b) show ojn dependence of the real 
and imaginary parts of the cumulant M{v,ujn), respectively, 
where r = (j, j) is the real-space vector connecting the sites 
i and j in the cluster [the inset of Fig. 14(b)]. We notice that 
r = (0, 0), (1,0), and (1,1) components are much larger than 
the other components at longer distances. This can be more 
clearly seen in Fig. 14(c), where the cumulant at the lowest 
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FIG. 15. (color online). Spectral functions at knode, calculated with 
iVc = 8 CDMFT+CTQMC fort' = -0.2t, U = 8t, and n = 0.95 at 
T = 0.05t (solid curve) and O.li (dashed curve), respectively. The 
arrow indicates the gap discussed in the text. Inset: The momenta 
(circles) specifying the superlattice consisting of the 8-site clusters. 



Matsubara frequency, ujq, is plotted against the Euclidean dis- 
tance, \Ji'^ + j^. We find that the cumulant is well localized 
within the 2x2 cluster, which justifies the M periodization. 

Although the cumulant at T = might be more extended 
than that at T — O.li, we do not have a tractable way to 
examine it for a cluster larger than 2x2. It is worthwhile, 
however, to note that the finite-temperature results obtained 
by our CDMFTh-CTQMC calculation are consistent with the 
CDMFTh-ED results with finite ry's, as discussed in Sec. IIIG. 
This fact implies that the M periodization with 2x2 cluster 



still remains a good approximation even at T = 0. 



APPENDIX C: NODAL SPECTRA OBTAINED BY 8-SITE 
CDMFT+CTQMC 

To demonstrate the fully-opened pseudogap in a larger 
cluster, we implement an 8-site cluster calculation with 
CDMFT+CTQMC. We calculate the spectral function 
A(k, oj) at k = knode = ('"■/S, 7r/2), where the d- and 
s-wave (fully-opened) pseudogaps are most distinguishable. 
The spectral function is obtained from G(k, iaj„) through the 
maximum-entropy method. 

Figure 15 shows A(k„ode, w) for t' = —0.2t, U — 8t, and 
n = 0.95. The small but clear gap just above the Fermi level 
at T = 0.05t proves the fully-opened pseudogap. The gap 
vanishes at a higher temperature T = O.lt. This is why the 
previous DCA study^^ at a high temperature T = 0.1 2i did 
not observe the gap. Note that, a very recent study with 8-site 
DCA^^ at T = O.Obt also detected a wispy reduction of the 
spectral weight at knode at low doping though the signal was 
so weak that it was not explicitly analyzed in their paper. The 
weakness of the signal would be due to the flat momentum 
dependence of the self-energy in each momentum patch in the 
DCA. Although the suppression of the spectrum is more sub- 
tle at higher doping in their result, it is reasonable because the 
pseudogap temperature decreases with doping. 

The gap formation at the nodal point cannot be reproduced 
by the assumption of the d-wave pseudogap. We have con- 
firmed this by employing the assumed form of the Green func- 
tion employed in Ref. 33. 
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